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Consistent Probabilistic Description of the Neutral Kaon System: Novel Observable 

Effects 
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Abstract 

We describe the neutral Kaon system as an open Lindblad-type quantum mechanical system due 
to Kaon decays. Including CP violation in the mixing, the formalism allows for total probability 
conservation, in contrast to the standard Weisskopf-Wigner Approach (WWA) employing a non- 
hermitian Hamiltonian to integrate out the decay products. This approach, in terms of density 
matrices for initial and final states, provides a consistent probabilistic treatment. Even for an initial 
pure state, the time evolution generates a mixed state, not described by a projector. To restore 
Unitarity, we include the dominant decay channel to two pions, so that one of the Kaon states 
with definite lifetime becomes stable. The new dynamics for initial and final states modifies the 
observables to second order in the CP-violating parameter e, as compared to WWA. 

PACS numbers: 14.40.Df; 03.65.Yz 



Neutral Kaons is a fascinating physical system that, 
due to its peculiar and at the time paradoxical behaviour 
in many respects, has lead to important discoveries, 
thereby triggering an enormous interest for its study. It 
is the first physical system where CP violation has been 
observed in the two-pion K° — >• 2 ir decay channel [lj , 
with the relevant experimental studies continuing up to 
date and extended to entangled neutral Kaon states 
in meson <fi factories Moreover, neutral Kaons have 
also been used as a probe of fundamental symmetries, 
such as CPT invariance [|| , and deviations from the stan- 
dard quantum mechanical behaviour. The latter may be 
induced by quantum gravity fluctuations appearing as a 
"decoherening" environment, leading to an open system 
(Lindblad-type Q) formulation [§-[l2j]. 

The standard description of the neutral-Kaon system 
follows the Weisskopf-Wigner Approach (WWA) [13[ for 
unstable particles using the non-Hermitian Hamiltonian 



H = M 



(1) 



other system with a Hamiltonian which is an abnormal 
operator. 

The WWA, restricted to the description of the dynam- 
ics of the initial two-state system, is furthermore unsat- 
isfactory from the point of view of securing Unitarity, 
i.e., conservation of probability. The latter property has 
to be imposed externally to the formalism, via the Bell- 
Steinberger Unitarity relation [ljj]. For this purpose, the 
inclusion of the tttt final state in the neutral-Kaon decays 
constitutes an excellent approximation towards restora- 
tion of Unitarity. 

The lack of a proper probabilistic interpretation of 
the neutral Kaon system has been addressed previ- 
ously [HI, [3] by distinguishing the ket and bra Hilbert 
spaces, so that the lack of orthogonality of the Kl, Kg 
states is bypassed. The treatment is, however, limited to 
considering the dynamics of pure initial states, without 
the inclusion of the final ones nor the evolution to mixed 
states. 
In 



17 



However, the simultaneous presence of CP violation in 
the mass matrix M and a difference of lifetimes in the 
antihermitian matrix iY /2 leads to a quantum incompat- 
ibility between M and T, [M,T\ ^ 0, i.e., one cannot 
define states of definite mass and lifetime simultaneously, 
because H is not a normal operator. The eigenstates Kl 
and Ks, obtained by a non- unitary diagonalization of H, 
lack physical meaning and their non-orthogonality pre- 
vents a consistent probabilistic treatment of this and any 



18] a suggestion has been made to view a de- 
caying quantum system as an open system interacting 
with an appropriate 11 environment" obtained by enlarg- 
ing the original Hilbert space by states representing the 
decay products. The time evolution of such a system 
can be described by an effective hermitian Hamiltonian, 
essentially M above, and an additional dissipative term 
of Lindblad form {dissipator) Q. As shown in the 
non-hermitian part of the Hamiltonian in the WWA, as- 
sociated with the particle decay width operator T, can 
be incorporated into the dissipator of the enlarged space 
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via a specific Lindblad operator B. 

To understand in simple terms the logic behind this 
open-quantum-system formalism for decaying systems, 
we first concentrate our attention on the evolution equa- 
tion for the initial density matrix, p — \^>)(^\ : 



p = —iH p + ipW = —i 



M,p 



(2) 



where p denotes time derivative. This equation can be 
formally obtained from the Schrodinger equation for the 
state vector |\&) with the non-Hermitian Hamiltonian %. 
In eq. ([2]), L is viewed as a single quantum mechanical 
operator. The non-Hermiticity of ~H leads to the anti- 
commutator term in the right-hand-side of this evolu- 
tion equation. As a consequence, the description of the 
system in terms of pure states, for which Trp 2 = Trp, 
breaks down at time t > 0. This can be readily shown 
by calculating the rate of the Von-Neumann entropy 
S = — Tr(plnp) using the evolution @, where only the 
anticommutator part contributes. Not only the system 
is evolving with time to mixed states, but Trp(i) ^ 0, so 
that the restriction to the initial Hilbert space, ignoring 
the decay products, leads to violation of Unitarity with 
non-conservation of Probability. 

To restore Unitarity we must include the final states 
when taking the trace Trp through a mapping from the 
initial Hilbert space to the final one (decay products) : 
Hi — > Hf. This mapping is implemented [18| by the 
transition operator B, which is related to T via: 



B^B = f 



(3) 



If {fk} denotes an orthonormal basis in Hf, and {ifij} 
denotes the corresponding orthonormal basis in Hi (or- 
thogonal to {fk} ), then one may write: 



(4) 



k=l 3 = 1 



where df = dimHf, and di = dimHi. The width op- 
erator r is thus a positive definite self-adjoint operator 
with non-negative eigenvalues. The latter can include 
possible zero eigenvalues, corresponding to stable states. 
This is to be contrasted to the corresponding expression 
given in [l8[ and will have important consequences for 
the neutral-Kaon system. 

The operator B can be considered as a sort of "en- 
vironment" operator from the point of view of the ini- 
tial state Hilbert space, and the evolution @ can be 
replaced now by an appropriate Lindblad evolution Q, 
with p spanning the combined initial (Hi) and final (Hf) 
Hilbert spaces, H to t = Hi@Hf. The Lindblad evolu- 
tion can be understood as follows: in view of <j3j> , the 
simple commutator structure of ([2]) in the conventional 
WWA [Toj j will now be replaced by an appropriate quan- 
tum ordering of the constituent operators B, B^ and p 



in such a way that the time evolution has the following 
properties 0] : (i) preserves the complete positivity of the 
density matrix operators at any time, i.e., the fact that 
their eigenvalues are positive or zero, so that the concept 
of probabilities associated with the eigenvalues of these 
operators makes sense, (ii) ensures the conservation of 
the total probability through Trp = 1 , including the final 
states (decay products) and (iii) implies increase of the 
entropy (of quantum states). 

Our density matrix p in the total Hilbert space H to t = 
Hiff)Hf is: p = I ^" ' ), where Hermiticity of p 

V Pfi Pff J 

is fulfilled in blocks and the subindices ii' (//') run over 
the initial (final) states. We have for the dimension of 
the relevant Hilbert spaces: 

dim Hi < oo , and dimHf >r = dimHi — no , (5) 

with no the degeneracy of the eigenvalue zero of the width 
operator. The evolution equations for the density matrix 
p in the H tot Hilbert space are then described by the 
Lindblad form Q: 



p = -i[H,p]-- (B^Bp + pB^B - 2BpB^) , (6) 



with 



U = H 1 



M 




B 




B 



(7) 



in total Hilbert space. The new formulation of the time 
evolution on the enlarged space has a hermitian Hamil- 
tonian and is probability conserving. The complete posi- 
tivity, that is guaranteed by construction in the Lindblad 
formalism Q, ensures that this feature characterises the 
decaying quantum system, exactly as it happens in sys- 
tems with Hermitian Hamiltonians. 

We would like to discuss here the application of this 
Lindblad open-system formulation of particle decay to 
physically realistic systems, such as neutral Kaons, which 
are known to exhibit CP violation and non-zero width 
difference Ar ^ 0. Contrary to WWA and the dynamics 
given by eq. ([2]) , the open-system formalism is applicable 
in terms of the transition operator B ([7]) , irrespective of 
the commutativity of the composite Y operator with M. . 
In this respect, the Lindblad dynamics ^ for the decay 
is appropriate mainly for neutral Kaons. Other neutral 
mesons, such as B-B systems, are characterised by very 
small width differences between the physical eigenstates, 
practically Ar ~ 0, for which the non- Hermitian Hamil- 
tonian is a normal operator and the WWA is satisfactory. 

For this discussion, we focus our attention from now 
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on two-state unstable systems. We can write eq. (|6]) as 



Pii> Pit 

Pfi Pff 



^-i[M,pw] - -{r ,p u /} 
y ipfiM - \pfiT 




(8) 



One notes the following: (i) the new dynamical be- 
haviour of the final state coupled to the initial one, with 
effects which cannot be described by T only; (ii) the for- 
mally identical structure of the equation for the time evo- 
lution of the initial-state density matrix, which is uncou- 
pled to the final states, to that of eq. ©, as a result of the 
anticommutator {T,pai}; however, here the evolution of 
pure to mixed states, in the sense of Trp(i) 2 7^ Trp(i) at 
a time t > 0, is evident due to the effects of the Lindblad 
operator B; (iii) the uncoupled dynamical behaviour of 
Pif(t), so that it is consistent to take Pif(t) = 0, if there 
is no initial (t = 0) mixed component between initial and 
final states; this implies that the description of the time 
evolution and decay is expressed in terms of initial and 
final density matrices only. 

The reader should recall once more that in the Lind- 
blad approach, total probability conservation for the den- 
sity matrix, including the decay products, is guaran- 
teed by construction, i.e., Trp.^ + Trpjf = 1 for any 
t, so that Unitarity is implied by the simple relation: 
dTr^m = _dTr(p ff (t)) ^ 

dt dt 
ified explicitly for the solutions we obtain here for the 

case of the neutral Kaon system. 

For the neutral Kaon K° — K° system we incorporate 
properly CP violation and the dynamics of its dominant 
decay to two pion final states. We use the 1^1,2) basis 
for the initial Kaon states defined as: 

\K X ) = ±={\K tt ) - \K )) , \K a ) = ±={\K°) + \K }), (9) 

which, as we show below, is a convenient choice in which 
the width operator is diagonal. 

The existence of a dominant decay channel in the neu- 
tral Kaon system, implies, via eq. ([5]), that no — 1, which 
is correct, given that there is only one zero eigenvalue in 
the spectrum of T. It would correspond to the Kl state 
in the WWA within time scales of order of the Ks life- 
time. Ignoring CPT Violation and CP violation in the 
decay, the choice of a real B leads to the result that the 
K1.2 states are the ones with definite lifetimes, so that 
the width operator in the \K\ t2 ) basis is given by the 
following 2x2 diagonal matrix with eigenvalues and 

r s : 



Fwwa 



i?e(r 12 ) 

r + i? e (r 12 ; 



= r s 




1 



(10) 



In this case the Lindblad operator ©, related to 
F via ©, is given by the following row matrix: 



B = ^/Ts~(Q , 1). In the \Ki, 2 ) basis ©, the mass M 
matrix, which will plav the role of the Hermitian Hamil- 
tonian, is written as [10(: 



M = 



M - Re(M 12 ) -Um{M 12 ) 
Um{M 12 ) M + Re(M 12 ) 



(11) 



ignoring again possible CPT- Violating effects. The CP 
violation parameter e is given by: 



|e|e- 



Im(M 12 ) 
+ iAm 



2Am , s 
tan<£= -— , (12) 



where Am = 2|M 12 | is the difference between the mass 
eigenvalues (A1.2 = M =F I-M12I) of the Kaon mass eigen- 
states. The latter are found to differ from the K\p 
states by terms of order of the CP violation parameter e: 
h) = \Ki) + 0{e)\K 2 ), \v 2 ) = \K 2 ) + 0(c)|Jfi). As this 
transformation between life-time and mass eigenstates is 
unitary (l9| . the orthogonality of the latter is guaranteed. 
This is in contrast to the WWA Kl and Ks states. 

To solve the evolution equations (jHJ, we shall follow 
the analysis in [1 01 ] . and use a perturbation method, by 
which we expand the density matrix elements at any time 
t in powers of the amplitude of the small CP-violation 
parameter |e| |12|: p u (t) = pf}(t) + p { j){t) + pf)(t) + 

. . . , p^) (t) oc |e| n , n = 0, 1, 2, . . . , where the indices /, J 
span the full Hilbert space of states {i, /}, including the 
decay products (final) states. 

In our analysis we shall restrict ourselves to second 
order in |e|, which matches the currently expected ex- 
perimental sensitivity. From equations ([5]). first one 
solves the evolution equation for the initial states pa'{t), 
i,i' = {1,2}, to order |e| 2 and then obtain pff(t), as- 
sociated with the / = (tt, it) decay channel, through 
Pff(t) = ^sP22{t)- The result, expressed in terms of 
the initial conditions for pu'(0), reads: 

P2 2 {t) = p 22 (0)e~ Fst 

-2|e||pi 2 (0)| [e~ rst cos(0 + 0i 2 ) - e~TT* cos(Ami - <\> - <p 12 )] 



[pn(0) + e rs * (pn(0) + p 22 (0)(2cos(2<rt + r s i)) 

(pn(0) cos(Amt) + ,922(0) cos(Amt - 20))] (13) 



-2e" 



with 4> 12 = Argp 12 (0). 

We can use this result in order to calculate various ob- 
servables of the Kaon system, in the above approximation 
of non-decaying K\ state. We can build useful observ- 
ables associated to the decay to irn or semileptonic decays 
ttIis. For a complete set of observables we refer the reader 
to a forthcoming publication [19]. For our purposes here 
we shall concentrate on two specific observables, namely 
the decay rates R(K° — > tttt) and R(K — > 7T7t), which 
will be compared to the corresponding observables in 
the WWA formalism, under the approximation Tl = 0, 
which is a good one for the range of times we consider. To 
be more specific, we will construct separately the sum of 
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(2(0 {io~ 5 ) 




FIG. 1: The quantity Q(t) = ^ £x 1 versus time (in 

-i j -v \ / ^ K|WWA 

units of the life-time rs), where T.R is given in eq. (|15[1 . 
The suffix WWA indicates the corresponding quantities in 
the WWA. 

rates, sensitive to even powers of |e|, and their difference 
(or CP-violating asymmetry), sensitive to odd powers of 

To this end we need the initially pure K°a.nd instates, 
prepared experimentally, in the 2) basis. These are 
described at t — in the total Hilbert space H tot by the 
density matrices: 

1 / 1 1 °\ 1 / 1 1 °\ 

Pk» = 5 1 1 , Pit o = - -1 1 . 

/ \000/ z \ / 

(14) 

To order \e\ we obtain: 

AR = R(K° -> tttt) - 7?(K -» tttt) = 
-2|e| r s e~ trs (cos - e"^' cos(Ami - 0)) 

Si? = i?(A"° -> tttt) + i?(A° -> tttt) = 
r s e" trs [l + \e\ 2 (1 + e trs + tT s + 

2 cos(20) - 4e~* cos( Ami - 0) cos </>) . (15) 

Comparing with the WWA formalism we observe that 
the quantities (|15l) differ from the corresponding expres- 
sions in WWA 2] by terms of order |e 2 |. The CP- 
violating decay-rate difference, AR, on the other hand, 
is of order |e|, agrees with the WWA result. In FIG. Q] 
we plot the difference from unity of the ratio of the sum 
of the rates for our result (fT5|) over that of the WWA as 
a function of time. We observe that differences of order 
10~ 5 occur already at times of order 5rg, which makes 
the experimental detection of the corrections challenging. 

For those readers concerned by the possible competi- 
tion of these novel effects of order |e| 2 with those coming 
from direct CP violation or order e' we point out the 
following: (i) The difference Q(t) is time dependent, (ii) 
what we call " (tttv) " in this work denotes the combination 
of rates ^ [2(7r + 7r~) + (7r° 7r )] , in which the contributions 
linear in e' cancel out. 

To conclude, we have presented a description of the de- 
caying neutral Kaon system as an open Lindblad system 



involving evolution of pure to mixed states. It satisfies 
all the physical requirements of a probabilistic quantum 
mechanical interpretation and guarantees Unitarity, pro- 
vided that the width operator in the dynamics of the 
initial states is a composite operator expressed in terms 
of the transition operator B between initial and final 
Hilbert spaces. As a consequence, the time dependence 
of the observable decay rates is modified, with respect to 
Weisskopf-Wigner Approach, to quadratic order in the 
magnitude of the CP violating parameter |e| 2 . In our 
view, the measurement of these novel observable effects 
constitute a new experimental challenge in the field. 
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